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1.0  INTRODUCTION 


Results  of  experimental  investigations  of  stability  characteristics  of  hypersonic  laminar  boundary 
layers  on  planar  and  conical  (sharp  and  blunt-nosed)  models  (Ref.l)  indicated  that  stability 
theory  can  form  a  useful  basis  for  describing  and  interpreting  wind  tunnel  data.  The  results  of 
(Ref.l)  indicate  that  the  measured  boundary  layer  stability  characteristics  of  simple 
configurations  were  described  well  in  a  general  sense  by  boundary  layer  stability  theory. 
Measurements  showed  well-defined  disturbances  in  the  laminar  boundary  layer  which  amplified 
until  breaking  down  into  turbulence.  Although  linear  stability  theory  was  only  partially 
successful  in  predicting  the  amplification  rate  of  these  disturbances,  their  frequency  range  and 
phase  velocity  was  predicted  well  by  linear  stability  theory.  Although  the  transition  Reynolds 
numbers  measured  in  conventional  wind  tunnels  are  generally  lower  than  in  free  flight 
(Ref.2),  the  results  of  Stetson  and  Kimmel  indicate  that  the  stability  process  in  a  hypersonic 
facility  like  AEDC  VKF  Tunnel  B  is  similar  to  the  presumed  free-flight  stability  mechanism. 

The  ambient  wind  tunnel  noise  tends  to  be  at  lower  frequencies  than  the  dominant  boundary  layer 
instability,  and  the  interaction  between  them  is  very  limited  (Ref.l.,  Ref.3).  The  lower  wind 
tunnel  transition  Reynolds  number  thus  results,  presumably,  from  larger,  but  still  linear  initial 
disturbance  amplitudes. 

The  supposition  that  boundary  layer  stability  in  a  wind  tunnel  environment  can  be  described  by 
linear  stability  theory  for  simple  configurations  (at  least)  opens  the  possibility  of  using  linear 
stability  theory  for  transition  prediction  and  design  of  wind  tunnel  models.  An  additional 
implication  is  that  wind  tunnel  transition  results  may  be  related  to  flight  transition  results  by 
using  linear  stability  theory.  Transition  in  a  wind  tunnel  might,  in  effect,  be  “calibrated”  for  a 
limited  range  of  models  and  conditions. 

The  stability  of  the  laminar  boundary  layer  on  a  sharp,  7-degree  half  angle  cone  with  freestream 
Mach  numbers  of  6.0  and  8.0  was  numerically  investigated  using  the  linear  stability  code  e-Malik 
(Ref .4).  Four  cases  were  studied,  each  of  which  modeled  wind  tunnel  conditions  from  tests 
performed  at  the  AEDC  Von  Karman  Hypersonic  Tunnel  B.  The  cases  studied  were: 

Table  1  AEDC  Wind  Tunnel  Test  Conditions 


CASE# 

Moc 

To[R] 

Twal|/To 

1 

8.0 

1310.00 

Adiabatic 

2 

8.0 

1310.00 

0.41 

3 

6.0 

1049.67 

0.63 

4 

6.0 

1310.00 

0.41 

Hot-wire  anemometer  measurements  were  used  to  determine  the  most  amplified  frequency  at 
transition.  The  comparison  of  measured  and  computed  most-amplified  frequency  provides  an 
indication  of  how  well  experimental  stability  processes  were  described  by  linear  stability  theory. 
Linear  stability  theory  was  used  to  compute  an  N-factor  to  predict  transition  for  each  of  the  cases. 
The  N-factor,  described  in  Section  4,  is  a  correlating  parameter  which,  when  measured  for  one 
case,  may  be  used  to  predict  transition  for  other  cases.  Case  1  provided  a  baseline  N-factor 


1 


which  was  then  used  to  predict  transition  for  the  other  cases.  Cases  2,  3,  and  4  provided  an 
assessment  of  wall  cooling  effects,  and  cases  2  and  4  provided  an  assessment  of  Mach  number 
effects. 
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2.0  Experimental  Configuration 


The  experimental  cases  are  documented  in  Refs.  (5,  6, 7)  and  summarized  in  Table  1. 

Boundary  layer  stability  data  were  obtained  using  hot-wire  anemometers  for  cases  1, 2,  and  4.  In 
addition,  heat  transfer  measurements  were  used  to  obtain  the  mean  transition  location  for  the 
cold- wall  cases  (2,  3  and  4).  Transition  is  determined  from  heat  transfer  measurements  by  noting 
where  heat  transfer  rates  first  rise  above  laminar  values.  Transition  was  measured  using  hot-wire 
anemometry  for  the  adiabatic  wall  case  1  by  noting  where  disturbance  levels  had  maximized  and 
begun  to  decay.  Previous  comparisons  of  hotwire  to  heat  transfer  measurements  indicate  a  good 
correlation  between  transition  locations  measured  with  these  two  techniques  (Ref.  8).  Figure  A1 
shows  the  experimental  setup  (hotwire  probe  and  sharp  cone)  at  AEDC  Tunnel  B. 

In  supersonic  flow  the  hot-wire  responds  to  fluctuations  of  mass  flux  and  total  temperature  (the 
sensing  variables)  whereas  the  fluctuations  which  characterize  an  unsteady  flowfield  are  those  of 
vorticity,  entropy  and  sound  (the  independent  mode  variables).  Fluctuations  in  the  sensing  and 
mode  variables  are  inferred  from  hot-wire  a-c  voltage  response  (rms.)  measurements  using  modal 
analysis.  For  detailed  information  on  modal  analysis  refer  to  (Ref.  1,12). 

At  each  model  station  the  hot-wire  probe  was  traversed  across  the  boundary  layer  to  locate  the 
point  of  maximum  rms  voltage  output  of  the  wire  (See  figure  A2).  The  quantitative  hot-wire 
data  were  acquired  at  these  maximum  energy  points  at  each  of  the  profile  survey  stations.  The 
recorded  signals  were  processed  using  the  spectrum  analyzer  followed  by  application  of 
computer  routines  to  eliminate  the  electronic  noise  and  to  restore  the  spectrum  of  the  wire 
response  data.  The  results  were  a  “noise-free”,  restored  amplitude-squared  spectral  density 
distribution  A^(f)  over  an  effective  frequency  range  from  zero  to  600  [kHz]  (See  figures  A3,  A4). 
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3.0  COMPUTER  CODE  e  -  Malik 


This  computer  program  (Ref.4)  solves  the  Stability  Equations  using  Spatial  Theory  (amplitude 
of  the  disturbance  changes  with  respect  to  space  not  time).  The  disturbance  wave  number  ( P )  in 
the  spanwise  direction  and  the  frequency  ( 0) )  are  prescribed  as  real  quantities.  The  numerical 
problem  is  to  obtain  the  eigenvalue  (oc ).  This  is  a  complex  quantity  indicating  the  wave  number 
in  the  streamwise  direction.  This  program  allows  the  user  to  perform  linear  stability  analysis  for 
2D  and  axisymmetric  flows  for  sharp  and  blunt  bodies.  It  also  allows  for  the  computation  of 
eigenfunctions  {w(y),v(y),/>(y),7’(y)}and  for  the  computation  of  N-Factor  curves  for  specified 

frequencies.  Section  3.0  describes  the  relation  between  eigenvalues,  eigenfunctions  and 
N-Factors. 
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4.0  SPATIAL  THEORY 

The  general  disturbance  in  a  boundary  layer  q'  (x,  y,  z,  t)  is  represented  as: 
q'  (x, y,  2,  t)  =  q(yy  exp[/(a5X  +  pgZ  - ttgcr)] 

where  the  non-dimensional  wave  number  in  the  freestream  (x)  direction  is: 

(X j  —  Otg^  -H  /(Xg|. 

the  non-dimensional  wave  number  in  the  spanwise  z-direction  is: 

Pa  ~  Par  ■’■^Pa; 

the  disturbance  phase  velocity  is: 
c  =  c^+  ic- 

Substituting  the  expressions  for  a,  p,c  into  the  disturbance  function  we  obtain  (after  some 
algebraic  manipulations),  the  following: 

q'  (x,y,z,t)  =  ^(y>=exp[(ag^c,.  +ag,.c,)^r - (ag,.x  +  pg,.z)+  i{a^^x  +  Pg,z  +  -otg^c,])] 

The  eigenfunction  is  given  as, 

q{y)=kriy)+i<ii{y)\ 

^(y)  =  k(yK 

where  tj  =  tan"'  [  is  the  eigenfunction  phase  angle  ;  therefore 

KW)) 

<l{y)  =  |^(y)|*[cos(ti)  -t-  i  sin(Ti)] 
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Observing  the  amplitude  function,  for  a  given  time  t,  the  amplitude  of  the  disturbance  is  constant 
along  the  line  ag^x  +  jSg^z  =  const  ,  this  is  a  line  of  constant  phase  in  the  x-z  plane.  This  constant 


phase  line  advances  with  a  phase  velocity  c 


ph 


a. 


When  fig  =0 ,  the  direction  of  propagation  of  the  wave  is  in  the  freestream  or  x  direction.  If 
fig  ^0 ,  the  wave  propagation  direction  is  inclined  at  an  angle  y/  with  respect  to  the  freestream  x 
direction.  The  wave  propagation  direction  is  always  normal  to  the  constant  phase  line  in  the  x  -z 


plane.  The  angle  is  given  by;  yr  =  tan 


-1 


\OCgj 


.  Even  if  the  boundary  layer  basic  state  is  two 


dimensional,  oblique  disturbances  are  three  dimensional. 

For  the  spatial  theory  development  (Ref.9),  time  dependency  is  eliminated,  therefore  the 
amplitude  function  becomes: 

A  =  |<?(y)|*exp[-(ag,.x  +  fi^z)] 

For  a  two  dimensional  disturbance  (amplitude  only  changes  in  the  x  -  direction),  fig^  =0  therefore 


the  amplitude  function  becomes: 

A  =  |9(y)|*exp[-(asx)] 

The  change  in  amplitude  with  respect  to  x  variable  is  known  as  the  spatial  amplification  rate  and 
is  given  by  the  following  equation  : 

But  A  =  k(y)|*exp[-(a5.x)]  is  the  original  amplitude  function,  therefore, 

—  =  -a 5,  A ,  where  (-a^  )  is  known  as  the  amplification  factor. 


The  N-Factor  is  defined  as  follows  : 


— In  A  = 
dx 


A  dx  j 


N{f)  =  -jag,dx 


—  =  ;  were  An  is  the  initial  amplitude  of  the  disturbance  at  the  lower  neutral  bound  where 

the  disturbance  first  begins  to  amplify. 


The  N-Factor  is  a  measure  of  how  much  a  particular  disturbance  amplitude  grows  as  it  travels 
downstream  relative  to  its  initial  amplitude.  For  this  study,  it  was  computed  for  each  specific 
frequency  for  a  range  of  Reynolds  numbers.  The  computation  was  repeated  for  various 
frequencies  until  the  most  amplified  frequency  was  obtained  at  the  transition  Reynolds  number. 
Much  work  has  shown  that  an  N-factor  of  9  may  be  correlated  with  transition  in  a  variety  of 
free-flight  conditions  (Ref  1 1).  This  correlation  pre-supposes  many  conditions,  such  as  the 
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freestream  disturbance  levels  being  of  equal  magnitude  for  all  conditions,  and  that  non-linear 
effects  prior  to  transition  are  equal  for  all  cases,  and  may  be  lumped  in  with  the  linear  instability. 
Despite  these  assumptions,  the  e^  method  has  been  useful  in  rationally  taking  parameter 
variations  into  account  in  transition  prediction.  The  fundamental  step  in  using  e"^  in  non- 
freeflight  conditions  is  obtaining  a  value  of  N  which  correlates  transition.  Assuming  transition 
occurs  at  a  constant  terminal  disturbance  level  A,  the  N-value  will  be  lower  in  the  wind  tunnel 
environment  since  the  initial  Ao  value  will  be  higher.  In  this  study.  Case  1  was  used  to  determine 
the  N-factor  for  AEDC  Tunnel  B.  This  N-factor  was  then  used  to  predict  transition  on  cases  2,  3, 
and  4. 

Another  concept  that  needs  to  be  introduced  is  the  subject  of  normal  modes  of  oscillation  (1^'., 
2"'*,  ...,etc.).  Normal  modes  are  easier  to  visualize  in  the  case  of  solid  media.  The  normal  modes 
of  the  boundary  layer  are  exactly  analogous.  A  cantilever  beam  system  will  be  used  to 
demonstrate  this  concept.  In  general,  if  a  sinusoidal  excitation  is  applied  to  a  beam  at  a  resonant 
frequency  (this  system  has  more  than  one  resonant  frequency)  also  referred  to  as  a  normal  mode 
frequency,  and  then  removed,  the  instantaneous  deflection  of  the  beam  is  referred  to  as  the 
normal  mode  shape.  The  eigenfunctions  obtained  for  this  problem  give  the  normal  mode  shapes. 
Like  the  cantilever  beam,  the  boundary  layer  is  an  oscillator  with  resonance’s.  Instability  waves 
correspond  to  normal  modes  of  the  boundary  layer  and  although  they  are  more  difficult  to 
identify  from  calculations  and  measurements  than  in  solid  mechanical  systems,  they  consist  of 
well  defined  disturbances  with  particular  mode  shapes,  or  eigenfunctions  and  specific 
combinations  of  phase  velocity  and  frequency. 

When  the  mean  flow  relative  to  the  disturbance  phase  velocity  (Cs)  is  supersonic  over  some 
portion  of  the  boundary  layer  mean  velocity  profile,  there  are  an  infinite  number  of  wave 
numbers  (a)  (solutions  of  the  eigenvalue  problem)  for  a  single  phase  velocity.  These  additional 
disturbances  are  called  higher  modes  (Mack,  ref  9).  Freestream  Mach  number  and  wall  cooling 
dictates  which  mode  of  oscillation  will  dominate.  For  subsonic  flow  the  amplification  rates  are 
higher  for  first  mode  type  oscillations,  for  supersonic  flow  the  amplification  rates  are  higher  for 
second  mode  type  oscillations. 
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5.0  NUMERICAL  COMPUTATIONS 


Program  e-Malik  was  used  to  obtain  amplification  rates,  N-Factor  curves  and  eigenfunctions  for 
specified  frequencies.  The  input  conditions  to  the  code  matched  the  various  wind  tunnel 
conditions  used  in  tests  performed  at  AEDC  Tunnel  B.  These  results  were  later  used  to  compute 
perturbation  contours  of  {p' {x,y,t),T {x,y,t),m! {x,y,t)] . 

The  inputs  to  the  program  for  each  test  case  are  listed  in  Table  1. 


Table  2  Program  Input  Data 


CASE# 

To[R] 

Twall/To 

Rex  [ft]-' 

1 

8.0 

6.8 

1310.00 

Adiabatic 

1.0e06 

2 

8.0 

6.8 

1310.00 

0.41 

1.0e06 

3 

6.0 

5.2 

1049.67 

0.63  n 

1.0e06 

4 

6.0 

5.2 

1310.00 

0.41 

1.0e06 

For  the  computation  of  the  perturbation  contours,  program  e-Malik  was  used  to  compute  the 
eigenfunctions  at  the  transition  location.  This  location  was  obtained  during  the  wind  tunnel  tests 
by  using  heat  transfer  and  hot-wire  measurements  for  each  of  the  cases  studied.  The  perturbation 
contours  were  obtained  by  expanding  the  eigenfunctions  in  the  streamwise  direction  using  the 
following  equation, 

q'(R,y,t  =  0)  =  \q(y)\*e^ *cos  [ttg,  * +T|  ] 
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6.0  GENERAL  RESULTS 


Measured  (where  available)  and  computed  values  of  most-amplified  frequency  and  transition 
Reynolds  number  are  summarized  in  Table  3.  These  computations  are  based  on  2-D  disturbances 
only.  This  is  justified  since  2-D  second  mode  disturbances  have  the  highest  amplified  rates  at 
these  Mach  numbers  (Ref.  8). 


Table  3  Transition  frequencies  and  Reynolds  numbers 


Measured  Rex 

Predicted  Rex 

Measured  fmax 
IkHz] 

Predicted  fmax 
fkHzl 

adiabatic 

HEiSEHi 

90 

90 

2 

0.41 

110 

120 

3 

lEEMEBi 

125 

100 

4 

0.41 

■eesebi 

not  available 

145 

A  comparison  of  measured  transition  Reynolds  number  to  computed  N-factors  for  Case  1 
indicated  a  transition  N-factor  of  4.4  for  AEDC  VKF-B.  This  N-factor  value  was  then  used  to 
predict  transition  for  the  other  three  cases.  Figures  A5,  A6,  A7  and  A8  show  the  computed  N- 
factor  curves  for  various  frequencies. 

In  general,  the  maximum  amplified  frequency  (fmax)  at  transition  was  predicted  well  by  the  e^ 
method.  Stetson  and  Kimmel  (Ref.l)  have  shown  that  the  maximum  amplified  frequency  is 
equal  to  approximately  0.5Ue/  5  for  a  number  of  cases.  This  is  because  typically,  the  second 
mode  wavelength  is  approximately  2  5,  and  its  convection  velocity  is  nearly  equal  to  Ue  (Ref.9). 

The  predicted  transition  Reynolds  number  for  the  cold  wall  Mach  8  (case  2)  agrees  well  with 
the  measured  value.  This  indicates  that  the  basic  stability  mechanism  and  the  wall  temperature 
effects  on  it  are  described  well  by  linear  stability  theory.  Although  linear  stability  theory  was 
inaccurate  in  predicting  the  disturbance  amplification  rates  for  Case  1  (Ref.  10),  the  good 
agreement  between  the  predicted  and  measured  cold  wall  transition  Reynolds  numbers  in  Case  2 
indicates  that  any  errors  must  be  consistent  between  the  two  cases  (1  &  2).  Linear  stability 
theory  predicts  a  destabilizing  effect  of  reduced  wall  temperature  ratio,  with  increasing  second 
mode  amplification  rates.  Also,  the  maximum  amplified  frequency  increases  due  to  the  decrease 
in  boundary  layer  thickness  with  wall  cooling,  this  is  observed  when  comparing  Cases  1  and  2. 

The  transition  Reynolds  number  for  Case  3  was  overpredicted  by  the  e”^  method  by  18%, 
although  the  trend  of  lower  transition  Reynolds  number  compared  to  Case  1  is  accurately 
replicated.  Mack  (Ref.  8)  shows  that,  for  2D  mean  flows,  second  mode  amplification  rates  begin 
to  dominate  at  M>4.0.  Mack’s  description  is  for  general  2D  mean  flow.  He  also  shows  that  as 
freestream  Mach  number  increases,  the  boundary  layer  becomes  more  stable,(transition  Reynolds 
number  increases).  When  comparing  Case  3  to  Casel,  since  they  both  dominated  by  second 
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mode  waves,  a  lower  transition  Reynolds  number  is  to  be  expected  for  Case  3,  since  both  wall 
cooling  and  decreased  Mach  number  destabilize  the  second  mode  boundary  layer. 

Case  4  is  identical  to  Case  3,  except  that  the  stagnation  temperature  is  increased  and  the  wall 
cooling  is  increased,  to  produce  a  ratio  of  (Tw  /  To)  of  0.41.  The  transition  Reynolds  number  was 
underpredicted  by  20%.  The  e^  method  predicts  a  transition  Reynolds  number  of  2.9x10^,  about 
73%  of  the  predicted  Case  3  transition  Reynolds  number.  Several  competing  effects  combine  to 
create  this  lower  predicted  transition  Reynolds  number.  First,  the  increased  wall  cooling 
destabilizes  the  second  mode,  but  stabilizes  the  first  mode.  Since  second  mode  disturbances  first 
begin  to  amplify  as  first  mode,  this  would  have  an  effect  on  the  maximum  amplitude  which  a 
disturbance  reaches.  Since  the  second  mode  growth  rates  are  much  higher  than  the  first, 
however,  the  effect  of  wall  cooling  on  the  second  mode  would  be  expected  to  dominate  the  effect 
of  wall  cooling  on  the  first  mode.  Second,  increased  stagnation  temperature  stabilizes  the 
boundary  layer  (Ref.  12).  The  measured  transition  Reynolds  number,  however,  is  almost  identical 
to  Case  3.  A  measured  maximum  amplified  frequency  is  not  available  for  this  case,  but  the 
computed  maximum  amplified  frequency  was  145  kHz.  This  is  significantly  higher  than  any  of 
the  other  cases  examined.  This  high  frequency  is  due  to  the  combination  of  a  thinner  boundary 
layer,  due  to  the  increased  wall  cooling,  and  a  higher  edge  velocity,  due  to  the  increased 
stagnation  temperature.  This  high  frequency  indicates  that  perhaps  a  third  effect  is  coming  into 
play  which  is  not  accounted  for  in  this  application  of  the  e*^  method.  Previous  measurements  of 
the  ambient  wind  tunnel  noise  in  this  facility  (Ref  1)  indicated  that  the  amplitude  of  freestream 
disturbances  drops  off  with  increasing  frequency,  roughly  proportional  to  1  /  f  When  the 
predicted  maximum  amplified  frequency  becomes  sufficiently  different  from  the  baseline  case 
used  to  determine  an  N-factor,  the  initial  Ao  used  to  describe  the  N-factor  changes.  In  this  case, 
since  the  predicted  most-amplified  frequency  is  much  higher  than  the  most-amplified  frequency 
of  the  baseline  case,  the  initial  Ao  amplitude  is  lower,  leading  to  a  higher  N-factor  at  transition, 
assuming  that  A(transition)  is  identical  for  all  cases.  It  would,  perhaps,  be  appropriate  to  take 
this  into  account  in  future  transition  predictions. 

Figure  A1  shows  the  AEDC  Tunnel  B  Experimental  Setup.  7-degree  sharp  angle  cone  and 
single  hot-wire  assembly. 

Figure  A2  shows  the  hot-wire  anemometer  peak  fluctuation  energy  signal  located  at  a  height  of 
85-90%  of  the  boundary  layer  thickness  (6)  for  M_=8.0  flow.  Hot-wire  sensitivity  coefficients  are 
a  function  of  the  mean  flow.  If  the  hot-wire  is  always  positioned  at  a  location  where  the  mean 
flow  is  the  same,  then  the  sensitivity  coefficients  do  not  change.  For  similar  boundary-layer 
profiles  (e.g.  sharp  cone)  if  the  wire  is  placed  at  the  same  y/5  position  in  the  boundary  layer,  as 
the  X  station  along  the  model  is  varied,  disturbance  amplification  rates  may  then  be  obtained 
directly  from  the  hot-wire  rms  signals. 

Figure  A3  shows  Fluctuation  spectra  at  the  location  of  peak  energy  in  the  boundary  layer  at 
approximately  0.95  for  adiabatic,  M^=8.0,  sharp  cone  flow.  Second  mode  disturbances  grow 
much  faster  than  first  mode  disturbances  and  rapidly  become  the  dominant  mode. 
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Figure  A4  contains  the  same  data  as  figure  A3  with  spectral  data  from  several  stations  overlaid 
to  better  illustrate  the  disturbance  frequencies.  The  first  and  second  mode  fluctuation  frequency 
bands  are  merged. 

Figures  A5  through  A8  show  the  computed  N-Factor  curves.  Case  #1  was  used  as  the  baseline 
for  obtaining  an  N-Factor  for  AEDC  Tunnel  B.  The  measured  transition  Reynolds  number  was 
4.7  X  10^,  the  corresponding  transition  N-Factor  value  obtained  from  figure  1  is  4.4.  As  stated 
before,  A/Ao  =  e^,  therefore,  AJAo  =  =  81.45.  The  amplitude  of  an  initial  disturbance  has 

increased  81.45  times  its  initial  amplitude  (at  the  lower  neutral  stability  curve)  by  the  time  it 
reaches  transition.  This  value  of  N  =  4.4  was  used  to  obtain  the  transition  location  for  the 
remaining  cases  from  figures  A6,  A7  and  A8. 


Figures  A9  through  A20  show  the  computed  magnitude  of  the  eigenfunctions.  The  general 
disturbance  can  be  expressed  as  : 
q'=  Ae'®  =  A[cos(0) -h  1  sin(0)] 

A=|?(y)|c-‘" 

6  =  \a^x\<r-  phase ■  angle 

The  magnitude  of  the  eigenfunction  is  contained  in  the  amplitude  function  of  the  disturbance. 
All  quantities  (pressure,  temperature,mass)  are  described  by  this  equation.  For  a  specific  case, 
the  magnitudes  between  the  quantities  are  different  but  the  amplification  rate  is  the  same.  The 
shapes  of  the  eigenfunctions  indicate  where  the  strongest  signals  are  going  to  be  located  across 
the  boundary  layer.  For  hotwire  measurements  we  are  mostly  interested  in  the  massflux 
eigenfunction  (hotwires  are  sensitive  to  massflux  and  total  temperature).  Figures  A17  through 
A20,  show  a  maximum  away  from  the  wall  close  to  the  boundary  layer  edge.  This  is  why  the 
strongest  signals  are  measured  near  the  boundary  layer  edge. 


Figures  A21  through  A32  show  the  computed  eigenfunction  phase  angle.  The  eigenfunctions 
are  a  complex  quantities.  The  general  representation  is  : 

^(y)  =  9.  (>)  +  «■<?, (y) 

The  eigenfunction  phase  angle  is  given  by  : 


T7  =  tan" 


This  angle  is  contained  in  the  general  2D-disturbance  equation  given  as  : 
q' (x,  y,t  =  0)  =  \q{y)\e'“<''  cos(a,x H-  r]) 


•q  represents  a  phase  shift  in  the  general  cosine  function.  Pressure,  temperature  and  massflux 
phase  shifts  are  different  since  the  eigenfunctions  are  different.  The  phase  shift  would  be 
measured  in  a  wind  tunnel  experiment  by  placing  two  hot-wire  probes,  one  above  the  other  in  the 
boundary  layer,  and  measuring  the  time  delaying  the  signal  between  the  two  probes. 


Figures  A3  3  through  A44  show  the  2D  contours  of  the  real  part  of  the  eigenfunction.  These 
plots  show  how  the  fluctuations  begin  to  amplify  as  they  travel  downstream.  The  N-Factor 
computation  indicate  the  rate  of  amplification  of  the  disturbances.  The  equation  describing  these 
plots  is : 
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q'iR,y,t  =  0)  =  \q{y)\e^  cos{a^R  +  q) 

These  contours  are  generated  for  a  single  frequency.  The  frequency  chosen  for  each  case  was  the 
maximum  predicted  frequency  (Table  1).  The  N-curves  used  are  the  ones  computed  for  each 
maximum  predicted  frequency,  a  different  curve  for  each  case  study. 

Figures  A45  through  A48  show  show  the  3D  representation  of  the  real  part  of  the  mass 
eigenfunction.  The  equation  used  to  generate  them  is  the  same  as  the  one  used  to  generate  figure 
A37  through  A40.  The  most  amplified  frequency  for  each  case  is  represented. 

Figures  A49  through  A52  show  the  real  part  of  the  mass  eigenfunction  at  a  constant  y  location. 
The  equation  used  to  generate  them  is  the  same  as  the  one  used  to  generate  figure  A41  -  A44. 

The  y  location  chosen  is  where  the  magnitude  of  the  eigenfunction  curve  shows  a  maximum.  See 
figures  A 17  through  A20. 
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7.0  CONCLUSIONS 


The  comparison  study  between  cases  1  and  2  and  between  cases  3  and  4  showed  that  the 
stability  parameters  are  very  sensitive  to  the  ratio  of  (Tw  /  To).  The  smaller  ratio  of  (Tw  /  To) 
resulted  in  a  decrease  in  the  predicted  transition  Reynolds  number.  This  indicates  that  wall 
cooling  destabilizes  second  mode  waves  by  increasing  their  amplification  rate  and  causing 
transition  to  occur  more  rapidly  (Cases  3  and  4). 

Freestream  Mach  number  dictates  what  mode  (first,  second)  will  be  amplified  the  most,  once 
the  dominant  mode  is  established,  amplification  rates  of  particular  frequencies  become  sensitive 
to  thermal  conditions,  specifically  to  the  ratio  of  (Tw  /  T©). 

Linear  stability  theory  properly  described  the  processes  leading  towards  transition  at  AEDC 
Tunnel  B  by  matching  ( Case  1)  or  coming  close  to  matching  the  predicted  maximum  amplified 
frequency  with  the  measured  frequencies.  This  indicates  that  AEDC  Tunnel  B  can  be  used  to 
conduct  stability  tests  and  produce  reliable  test  data  with  the  use  of  transition  N-Factor 
correlation’s  for  given  geometry’s  and  test  conditions. 
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APPENDIX  A 


Figure  A 1  Experimental  Setup  AEDC  Tunnel  B 
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Figure  A  4  Fluctuation  spectra  overlaid 
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Figure  A  5  N-Factor  Curves,  M=8.0,  Adiabatic 
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N~Factor  vs.  SqrtfRey] 

Me=5.2 ;  To=l049.67[R] ;  Tw/ro=0.63 
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Figure  A  7  N-Factor  Curves,  M=6.0,  Tw/To=0.63 
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Figure  A  8  N-Factor  Curves,  M=6.0,  Tw/To=0.41 
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Figure  A  9  Mag.  Pressure  eigenfunction,  M=8.0,Adiabatic 


Figure  A  10  Mag.  Pressure  eigenfunction,  M=8.0,  Tw/To=0.41 
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Figure  All  Mag.  Pressure  eigenfunction,  M=6.0,  Tw/To=0.63 


Figure  A 12  Mag.  Pressure  eigenfunction,  M=6.0,  Tw/To=0.41 
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Figure  A 13  Temperature  eigenfunction,  M=8.0,  Adiabatic 


Figure  A  14  Temperature  eigenfunction,  M=8.0,  Tw/To=0.41 
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Surface  Normal  (y) 


Figure  A  15  Temperature  eigenfunction,  M=6.0,Tw/To=0.63 


Figure  A 16  Temperature  eigenfunction  M=6.0,  Tw/To=0.41 


Figure  A 17.  Mass  eigenfunction,  M=8.0,  Adiabatic 


Figure  A  18  Mass  eigenfunction,  M=8.0,  Tw/To=0.41 
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Figure  A 19  Mass  eigenfunction,  M=6.0,  Tw/To=0.63 


Magnitude  Mass  Eigen  Function 
T  =1310.0[R] ;  M=6.0  Cold  Wall 
T.«/T,  =  0.41 


Figure  A  20  Mass  eigenfunction,  M=6.0,  Tw/To=0.41 
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Figure  A  22  Pressure  function  phase  angle,  M=8.0,  Tw/To=0.41 
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Surface  Normal 


Figure  A  24  Pressure  function  phase  angle,  M=6.0,  Tw/To=0.41 
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Temperature  Eigen  Function  Phase  Angle 
T  =I310.0[R1 ;  M=8.0  Adiabatic 


Figure  A  25  Temperature  function  phase  angle,  M=8.0,  Adiabatic 
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Figure  A  26  Temperature  function  phase  angle,  M=8.0,  Tw/To  =  0.41 
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Surface  Normal  (y)  |  |  Surface  Normal 
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Figure  A  29  Mass  function  phase  angle,  M=8.0,  Adiabatic 


Figure  A  30  Mass  function  phase  angle,  M=8.0,  Tw/To  =  0.41 
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Real  Part  Pressure  Function  Re{p’(x,y,t=0)} 
T^=:13i0,0tR] ;  MrrS.O  Adiabadc 
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Figure  A  33  Real  part  pressure  function,  M=8.0,  Adiabatic 


Figure  A  34  Real  part  pressure  function,  M=8.0,  Tw/To  =  0.41 
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Surface  Normal  (y) 


Real  Part  Pressure  Function  Re[p’  (x,y,t=0)} 
T^=1049,ti7IR] ;  Cold  Wall 
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Figure  A  35  Real  part  pressure  function,  M=6.0,  Tw/To  =  0.63 
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Figure  A  36  Real  part  pressure  function,  M=6.0,  Tw/To=0.41 
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Figure  A  38  Real  part  temperature  function,  M=8.0,  Tw/To  =  0.41 
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Surface  Normally) 


Figure  A  39  Real  part  temperature  function,  M=6.0,  Tw/To  =  0.63 
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Figure  A  42  Real  part  mass  function,  M=8.0,  Tw/To  =  0.41 
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Surface  Normal  (y)  ^  I  Surface  Normalfy) 


Real  Part  Mass  Function  Rc{m’  (x,y,t3=0)  } 
T  =1049.67IR] ;  M=6.0  Cold  Wall 
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Figure  A  43  Real  part  mass  function,  M=6.0,  Tw/To  =  0.63 
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Real  Part  Mass  Funedon  Re{ni’(74y>t=0)} 
T^=1310.0[R] ;  M=6.0  Cold  Wall 
T^/T-0.41 


ReCm*<x,y,t»0)l 


207J224 

184,559 

J61.S94 

139.229 

116.564 

93.S992 

71.2343 

48.5695 

25.9046 

3.2396$ 

-19,4252 

-410901 

-64.755 

-874199 

-nac85 

-13175 

-155.415 

-178.079 

-200.744 

-223.409 


Figure  A  44  Real  part  mass  function,  M=6.0,  Tw/To  =  0.41 
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Real  Part  Mass  Function  Re{in^(x,y,t=0)} 
T^:rl310.0[R3 ;  M=::«.0  Adiabatic 
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Figure  A  45  Real  part  mass  function,  M=8.0,  Adiabatic 


Figure  A  46  Real  part  mass  function,  M=8.0,  Tw/To  =  0.41 
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Figure  A  47  Real  part  mass  function,  M=6.0,  Tw/To  =  0.63 


Figure  A  48  Real  part  mass  function,  M=6.0,  Tw/To  =  0.41 
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Figure  A  49  Real  part  mass  function,  M=8.0,  Adiabatic 


Figure  A  50  Real  part  mass  function,  M=6.0,  Tw/To  =  0.41 
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m*(x,yi=4*30,t=0)} 


Figure  A  52  Real  part  mass  function,  M=6.0,  Tw/To  =  0.41 
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